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O ■ Abstract 



It is shown that the previously noted extreme perturbative NNLO/NLO instabihty 
of the longitudinal structure function Fi{x, Q^) is a mere artefact of the commonly 
5^ , utilized 'standard' gluon distributions. In particular it is demonstrated that using 

the appropriate - dynamically generated - parton distributions at NLO and NNLO, 
Fl{x,Q'^) turns out to be perturbatively rather stable already for Q^ > O (2 — 3 
GeV2). 



A sensitive test of the reliability of perturbative QCD is provided by studying [1, 2, 3, 4] 
the perturbative stability of the longitudinal structure function Fi{x, Q"^) in the very small 
Bjorken-x region, x ^ 10~^, at the perturbatively relevant low values of Q^ ^ (9(2 — 3 
GeV^). For the perturbative-order independent rather flat toy model parton distributions 
in [1], assumed to be relevant at Q^ — 2 GeV^, it was shown that next-to-next-to-leading 
order (NNLO) effects are quite dramatic aX x ^ 10^^ (cf. Fig. 4 of [1]). To some extent 
such an enhancement is related to the fact, as will be discussed in more detail below, 

that the third-order a^ contributions to the longitudinal coefficient functions behave like 

(s) 
xc\ ~ — Inx at small x, as compared to the small and constant coefficient functions at 

LO and NLO, respectively. It was furthermore pointed out, however, that at higher values 

of Q^, say Q^ ~ 30 GeV^, where the parton distributions are expected to be steeper in 

the small-x region (cf. eq. (13) of [1]), the NNLO effects are reduced considerably. It 

is well known that dynamically generated parton distributions [5] are quite steep in the 

very small-z region already at rather low Q^, and in fact steeper [6] than their common 

'standard' non-dynamical counterparts. Within this latter standard approach, a full 

NLO (2-loop) and NNLO (3-loop) analysis morevover confirmed [2, 3] the indications for 

a perturbative fixed-order instability observed in [1] in the low Q^ region; even additional 

resummations have been suggested [7] in order to remedy these instabilities. 

It is therefore interesting to study this issue concerning the perturbative stability of 
Fl{x,Q'^) in the low Q^ region, Q^ ^S 5 GeV^, within the framework of the dynamical 
parton model [5, 6]. For this purpose we repeat our previous [8] 'standard' evaluation 
of the NLO and NNLO distributions within the dynamical approach where the parton 
distributions at Q > 1 GeV are QCD radiatively generated from valenc^Wke (positive) 
input distributions at an optimally determined Qo = /^ < 1 GeV (where 'valence-like' 
refers to a/ > for a// input distributions xf{x, fi^) ~ x"'f{l — x)^f). This more restrictive 
ansatz, as compared to the standard approach, implies of course less uncertainties [6] 
concerning the behavior of the parton distributions in the small-z region aX Q > fi which 



is entirely due to QCD dynamics at a; <, 10 ^. The valence-like input distributions at 
Qo = fi < 1 are parametrized according to [8] 

xqv{x,Ql) = Ng^x"^- {1 - xf^- {1 + Cg^y/x + dg^x + egy-^) (1) 

xw{x,Ql) = Nyjx'''^ {1 - xf"^ {1 + c^y/x + d^x) (2) 

for the valence q^ = u^, d^ and sea w = q, g densities, and a vanishing strange sea at 
Q^ = Ql, s{x, Ql) = s{x, Ql) = 0. All further theoretical details relevant for analyzing F2 
at NLO and NNLO in the MS factorization scheme have been presented in [8], using again 
the QCD-PEGASUS program [9] for the NNLO Q^-evolutions, appropriately modified 
to account for the fixed Uf = 3 flavor number scheme with a running as{Q^). The 
heavy flavor (dominantly charm) contribution to F2 is taken as given by fixed-order NLO 
perturbation theory [10, 11] using rric = 1.3 GeV and nib = 4.2 GeV as implied by optimal 
fits [6] to recent deep inelastic c- and 6-production HERA data. Since a NNLO calculation 
of heavy quark production is not yet available, we have again used the same NLO 0(q;^) 
result. This is also common in the literature [12, 13, 14] and the error in the resulting 
parton distributions due to NNLO corrections to heavy quark production is expected [12] 
to be less than their experimental errors. Finally, we have used for our fit-analyses the 
same deep inelastic HERA-Hl, BCDMS and NMC data, with the appropriate cuts for F^'" 
as in [8] which amounts to a total of 740 data points. The required overall normalization 
factors of the data turned out to be 0.98 for HI and BCDMS, and 1.0 for NMC. We 
use here again solely deep inelastic scattering data since we are mainly interested in the 
small-x behavior of structure functions. The resulting parameters of the NLO and NNLO 
fits are summarized in Table 1. The dynamical gluon and sea distributions, evolved to 
some specific values of Q^ > Qo^ are at the NLO level very similar to the ones in [6] which 
were obtained from a global analysis including Tevatron Drell-Yan dimuon production 
and high-ET inclusive jet data as well. Furthermore, the dynamically generated gluon is 
steeper as a; ^ than the gluon distributions obtained from conventional 'standard' fits 
[6, 8] (based on some arbitrarily chosen input scale Qo > I GeV^, i.e. Qq — 2 GeV^). On 



the other hand, the dynamical sea distribution has a rather similar small-z dependence as 
the 'standard' ones [6, 8]; this is caused by the fact that the valence-like sea input in (2) 
vanishes very slowly as x ^ (corresponding to a small value of ttq, ttq ~ 0.07 according 
to Table 1) and thus is similarly increasing with decreasing x down to x ~ 0.01 as the 
sea input obtained by a 'standard' fit. Similar remarks hold when comparing dynamical 
and standard distributions at NNLO. At NNLO the gluon distribution xg is flatter as x 
decreases and, in general, falls below the NLO one in the small-x region, typically by 20 
- 30% at X ~ 10^^ and Q^ ^ 10 GeV^, whereas the NNLO sea distribution xq is about 
10 - 20% larger (steeper) than the NLO one. Furthermore it should be mentioned that 
the NLO q;s(M|) in Table 1 turns out be somewhat smaller in flts based solely on deep 
inelastic structure function data [8, 12, 15, 16, 17] as compared to those which take into 
account additional hard scattering data [2, 6, 18, 19] (for a recent summary, see [20]). 
At NNLO the resulting ^^(Ml) is generally slightly smaller [20] (c.f. Table 1) which is 
due to the fact that the higher the perturbative order the faster as(Q^) increases as Q^ 
decreases. 

Now we turn to the perturbative predictions for Fl{x, Q^) which can be written as 

2 

X-^Fl = CL,ns ®qns + - {C^q Qs + C^g ® Q) + X-^ Fl (3) 

y 

where ® in the n/ = 3 light quark flavor sector denotes the common convolution, Qns 
stands for the usual flavor non-singlet combination and g^ = 'Yliq=udsi'i + ^) ^^ ^^^ '^'-'^" 
responding flavor-singlet quark distribution. Again we use the NLO expression [10, 11] 
for F£ also in NNLO due to our ignorance of the 0{al) NNLO heavy quark corrections. 
The perturbative expansion of the coefficient functions can be written as 

0.,x) = 5^(^)"cg(x). (4) 

n=l ^ ^ 

In LO, Cj^ng = ^x, c\^ = 0, c^ = 24x(l — x) and the singlet-quark coefficient function 
is decomposed into the non-singlet and a 'pure singlet' contribution, c^^ = c^^g + c^^s- 



Sufficiently accurate simplified expressions for the exact [21, 22] NLO and [23] NNLO co- 
efficient functions c}^ ^ and c\^[, respectively, have been given in [1]. It has been futhermore 
noted in [1] that especially for Ci^g both the NLO and NNLO contributions are rather 
large over almost the entire x-range. Most striking, however, is the behavior of both CL,q 
and C^^g at very small values [1, 24] of x: the vanishingly small LO parts {xc\^{ ~ x^) are 
negligible as compared to the (negative) constant NLO 2-loop terms, which in turn are 

(3) 

completely overwhelmed by the positive NNLO 3-loop singular corrections xc^ ^ ~ — Inx. 
This latter singular contribution might be indicative for the perturbative instability at 
NNLO [1], as discussed at the beginning, but it should be kept in mind that a small-x 
information alone is insufficient for reliable estimates of the convolutions occurring in (3) 
when evaluating physical observables. 

For this latter reason we do not display the individual gluon and (sea) quark distribu- 
tions separately but instead our predictions for the convolutions of the individual light 
u, d, s quark (F^) and gluon (F^) contributions in (3) in Figs. 1 and 2, respectively, at 
two characteristic low values of Q^. (Note that Fl + F[ = Fl — F£ according to (3)). 
Although the perturbative instability of the subdominant quark contribution in Fig. 1 as 
obtained in a 'standard' ffi does not improve dramatically for the dynamical (sea) quark 
distributions, the instability disappears almost entirely for the dominant dynamical gluon 
contribution already at Q^ — 2 GeV^ as shown in Fig. 2. This implies that the dynamical 
predictions for the total Fl{x,Q'^) become perturbatively stable already at the relevant 
low values of Q^ ^ 0{2 — 3 GeV^) as evident from Fig. 3, in contrast to the 'standard' 
results in Fig. 4. In the latter case the stability has not been fully reached even aX Q^ = 5 
GeV^ where the NNLO result at x = 10"^ is more than 20% larger than the NLO one. 
A similar discrepancy prevails for the dynamical predictions in Fig. 3 at Q^ = 2 GeV^. 
This is, however, not too surprising since Q^ = 2 GeV^ represents somehow a borderline 
value for the leading twist-2 contribution to become dominant at small x values. This 
is further corroborated by the observation that the dynamical NLO twist-2 ffi slightly 



undershoots the HERA data for F2 at Q^ — 2 GeV^ in the small-z region (cf. Fig. 1 of 
[6]). The NLO/NNLO instabihties imphed by the standard fit resuhs obtained in [2, 3] 
at Q^ ~ 5 GeV^ are even more violent than the ones shown in Fig. 4. This is mainly due 
to the negative longitudinal cross section (negative Fl{x,Q'^)) encountered in [2, 3]. The 
perturbative stability in any scenario becomes in general better the larger Q'^, typically 
beyond 5 GeV^ [1, 2, 3], as shown in Figs. 3 and 4. This is due to the fact that the 
Q^-evolutions eventually force any parton distribution to become sufficiently steep in x. 

For completeness we finally compare in Fig. 5 our dynamical (leading twist) NNLO 
and NLO predictions for Fl{x, Q"^) with a representative selection of (partly preliminary) 
HERA-Hl data [25, 26]. Our results for Fl, being gluon dominated in the small-x region, 
are in full agreement with present measurements which is in contrast to expectations 
[2, 3] based on negative parton distributions and structure functions at small values of 
X. To illustrate the manifest positive definiteness of our dynamically generated structure 
functions at Q'^ > IJ-^ = 0.5 GeV^ we show Fl{x, Q'^) in Fig. 5 down to small values of Q^ 
although leading twist-2 predictions need not necessarily be confronted with data below, 
say, 2 GeV\ 

To summarize, we have shown that the extreme perturbative NNLO/NLO instability 
of the longitudinal structure function Fl at low Q^, noted in [2-4], is an artifact of the 
commonly utilized 'standard' gluon distributions rather than an indication of a genuine 
problem of perturbative QCD. In fact we have demonstrated that these extreme instabil- 
ities are reduced considerably already at Q^ = 2 — 3 GeV^ when utilizing the appropriate, 
dynamically generated, parton distributions at NLO and NNLO. These latter parton dis- 
tributions have been obtained from a NLO and NNLO analysis of Fl*'" data, employing 
the concepts of the dynamical parton model. It is gratifying to notice, once again, the 
advantage of the dynamical parton model approach to perturbative QCD. 
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Table 1: Parameter values of the NNLO and NLO QCD fits with the parameters of the input 
distributions referring to (1) and (2) at a common input scale Q^ = fj? = 0.5 GeV^ which turns 
out to be optimal at both perturbative orders. 
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Figure 1: The individual light {u,d,s) quark contribution F^ to the total F^ in (3) in 
the dynamical (dyn) and standard (std) parton approach at NNLO and NLO for two 
representative low values of Q^. The standard parton distributions utilized in the lower 
panel are taken from [8]. Our standard NLO results in the lower panel are similar for the 
CTEQ6 (anti)quark distributions [18]. Notice that, according to (3), F^ + F| = F^ — F£. 
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Figure 2: As in Fig. 1 but for the gluonic contribution Fl to F^ in (3) with Fl 
|x CL,g ® g. 



0.5 

0.4 



^ = 2 GeV^ 




TI, I I I 



i^L f-^' Q) 



I mil 1 — I I I iiiii 1 — I I I Mill 1 — I I I iiiii 1 — I I I II 



dyn 



NNLO 
NLO 



^ = 3 GeV^ 



1 1 iiiiii| — 1 1 iiiiii| — 1 1 iiiiii| — 1 1 iiiiii| — i-t+f 



0.4 
0.3 
0.2 



0.1 - 



(^ = 4 GeV^ 




0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 






e^ = 5 GeV^ 



1 1 iiiiii| — 1 1 iiiiii| — 1 1 iiiiii| — 1 1 iiiiii| — i-i+f 



(i = 10 GeV^ 



Q^ = 20 GeV^ 



111 I I I 



' ' ' ' ■^*'*-^""' ' I I ' ■^■'•♦-j-"" 



J I I I 



10"^ 10""^ 10"^ 10"^ 10"^ 10"^ 10""^ 10"^ 10"^ 10"^ 1 

X 



Figure 3: Dynamical parton model NNLO and NLO predictions for Fl{x, Q"^) in (3). 
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Figure 4: As in Fig. 3 but for the common standard parton distributions as taken, for 
example, from [8]. 
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Figure 5: Our dynamical NNLO and NLO predictions for Fl at a fixed value oiW = 276 
GeV. The (partly preliminary) HI data [25, 26] are at fixed W ~ 276 GeV. 



